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Abstract 

The spin-phase interference effects are studied analytically in resonant 
quantum tunneling of the Neel vector between degenerate excited levels in 
nanometer-scale single-domain antiferromagnets in the absence of an exter- 
nal magnetic field. We consider a model for mesoscopic antiferromagnets 
with uncompensated excess spins for the more general structure of magnetic 
anisotropy, such as biaxial, trigonal, tetragonal and hexagonal crystal symme- 
try. This study provides a nontrivial generalization of the Kramers degeneracy 
for double-well system to coherently spin tunneling at ground states as well 
as low-lying excited states in AFM system with m-fold rotational symmetry 
around the z axis. The energy level spectrum and the thermodynamic prop- 
erties of magnetic tunneling states are found to depend significantly on the 
parity of the excess spins at sufficiently low temperatures. Possible relevance 
to experiments is also discussed. 
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I. INTRODUCTION 



Can quantum mechanics be used to describe the behavior of macroscopic objects? This 
question has fascinated physicists for more than 70 years until the 1980s when it was pro- 
posed by Leggett et a/.§ that macroscopic objects could behave quantum mechanically, pro- 
vided that the dissipative interactions with the environment were small enough. In recent 
years, macroscopic quantum phenomena have been observed in various systems — for exam- 
ple, quantum tunneling of the phase in a Josephson junction, permanent current in small 
conductor rings, Bose-Einstein condensation in atomic vapours, and C 60 molecules. Recently, 
tunneling of the magnetization has been intensively studied theoretically and experimentally 
in nanoparticles and molecular clusters.! Two kinds of phenomena can be envisaged. First, 
macroscopic quantum tunneling (MQT), where the magnetization tunnels through an bar- 
rier from one metastable state to a stable one. Second, macroscopic quantum coherence 
(MQC), where the magnetization oscillates back and forth between the degenerate states. 
This oscillation should show up in the frequency-dependent magnetic noise and suscepti- 
bility spectraSi Besides its importance from a fundamental point of view, tunneling of 
the magnetization changes the properties of small magnets, with potential implications for 
the data-storage technology. It is also very important to the reliability of nanometer-scale 
magnetic units in memory devices and the designing of quantum computers.il 

One notable subject in spin quantum coherence is that the topological Wess-Zumino- 
Berry phasei can lead to remarkable spin-parity effects. It was found that the tunnel split- 
ting is suppressed to zero for half-integer total spins in biaxial ferromagnetic (FM) particles 
due to the destructive phase interference between topologically different tunneling pathsJl 
However, the phase interference is constructive for integer spins, and hence the splitting is 
nonzero.! While spin-parity effects are sometimes be related to Kramers degeneracy, they 
typically go beyond the Kramers theorem in a rather unexpected way.0 The auxiliary parti- 
cle method was proposed to study the model for a single large spin subject to the external 
and anisotropy fields, and to discuss the spin-parity effects.i Similar effect was found in 
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antiferromagnetic (AFM) particles, where only the integer excess spins can tunnel but not 
the half-integer ones.1'0 However, these spin-parity effects are intrinsically absent in the 
phenomena of MQT and MQC in the Josephson-j unction-based superconducting systems, 
which makes the tunneling phenomena in nanoscale magnets more important for under- 
standing the foundations of quantum mechanics. Theoretical results showed that MQC in 
AFM particles should show up at higher temperatures and higher frequencies than in FM 
particles of similar sizeJ^Hi! which makes AFM particles more interesting for experimen- 
tal purposes. Recently, topological phase interference effects were investigated extensively 
in FM and AFM particles in a magnetic field,@~i&lizl and in the systems with different 
symmetries.S'El 

One recent experiment^ was performed to measure the tunnel splittings in molecules 
Fes, and a clear oscillation of the splitting as a function of the field along the hard axis was 
observed, which is a direct evidence of the role of the topological spin phase (Berry phase) 
in the spin dynamics. Recent theoretical and experimental studies include the thermally 
activated resonant tunneling based on the exact diagonalization,ill the auxiliary particle 
method,! the discrete WKB method and a nonperturbation calculation,il the non-adiabatic 
Landau-Zener model ,@ and the calculation based on exact spin-coordinate correspondence.il 

The importance of the topological interference term of the Berry phase for the problem 
of spin tunneling and the associated spin-parity effects have been elucidated in Refs. 6, 7 
and 10. However, the theoretical studies on AFM systems§~@ have been focused on phase 
interference between two opposite winding ground-state tunneling paths in biaxial particles. 
The spin-phase interference between excited-level tunneling paths is not clearly shown for 
AFM particles. Moreover, the previous works on AFM spin tunnelinglHll have been confined 
to the system with biaxial symmetry, which has two energetically degenerate easy directions 
in the basal plane. The purpose of this paper is to study the quantum tunneling and 
spin-phase interference at excited states for AFM particles in the absence of an external 
magnetic field. By calculating the nonvacuum instantons, we obtain the analytical results 
for tunnel splittings at excited levels. To compare theory with experiment, we consider the 



AFM particles with the general structure of magnetocrystalline anisotropy, such as biaxial, 
trigonal, tetragonal, and hexagonal symmetry around z, which have two, three, four, and six 
energetically degenerate easy directions in the basal plane. For AFM particles with biaxial 
symmetry, the spin-phase interference effect can be studied by summing up the contributions 
of topologically different tunneling paths of clockwise and counterclockwise instantons. Mi 
However, for the system with complex trigonal, tetragonal or hexagonal symmetry, this 
procedure is hard to evaluate. In this paper, the spin tunneling problem is mapped onto 
a particle moving problem in one-dimensional periodic potential U (4>) by integrating out 
the momentum in the path integral, and the tunneling level spectrum of excited states 
is obtained by using the Bloch theorem. Our results show that the excited-level tunnel 
splittings depend significantly on the parity of the excess spins of AFM particles. The low- 
energy limits of the nonvacuum instanton and the tunnel splittings of excited levels agree 
well with the results of ground-state tunneling in biaxial AFM particles.0 The structure 
of tunneling level spectrum for the trigonal, tetragonal and hexagonal crystal symmetry is 
found to be much more complex than that for the biaxial crystal symmetrySHli The tunnel 
splitting for biaxial AFM particles is quenched to zero for half-integer excess spins due to 
the destructive interference of Berry phase.i&^ii While the tunnel splitting can be nonzero 
even if the excess spin is a half-integer for the trigonal, tetragonal, or hexagonal symmetry 
at zero magnetic field. 

This paper is structured in the following way. In Sec. II, we review briefly some basic 
ideas of quantum tunneling in AFM particles, and discuss the fundamentals concerning the 
computation of excited-level splittings in the double-well potential. In Sees. Ill we study 
the resonant quantum tunneling of Neel vector between degenerate excited states in AFM 
particles with the biaxial symmetry in detail, and present the results for trigonal, tetragonal 
and hexagonal symmetry in Sec. IV. The conclusions are presented in Sec. V. 
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II. PHYSICAL MODEL 



The system of interest is a single-domain AFM particle of about 5~10 nm in radius at a 
temperature well below its anisotropy gap. According to the two-sublattice model,0 there 
is a strong exchange energy mi • 1112/xi between two sublattices, where mi and rri2 are the 
magnetization vectors of the two sublattices with large, fixed and unequal magnitudes, and 
X± is the transverse susceptibility. Under the assumption that the exchange energy between 
two sublattices is much larger than the magnetocrystalline anisotropy energy, the Euclidean 
action for the AFM particle (neglecting dissipation with the environment) i 
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where 7 is the gyromagnetic ratio, a is the exchange constant (which is also referred to as the 
stiffness constant, or the Bloch wall coefficientl!) , and r = it is the imaginary-time variable. 
The E(6, (p) term includes the magnetocrystalline anisotropy and the Zeeman energies. The 
polar coordinate 9 and the azimuthal coordinate 0, which are the angular components of 
mi in the spherical coordinate system, can determine the direction of the Neel vector. 

As pointed out in Ref. 10, for a nanometer-scale single-domain AFM particle, the Neel 
vector may depend on the imaginary time but not on coordinates because the spatial deriva- 
tives in Eq. (1) are suppressed by the strong exchange interaction between two sublattices. 
So all the calculations performed in the present work are for the homogeneous Neel vector. 
Therefore, Eq. (1) reduces to 
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where V is the volume of the single-domain AFM nanoparticle. m = mi — m 2 = hjs/V, 
where s is the excess spin due to the noncompensation of two sublattices. Note that the first 



term in the Euclidean action is a total imaginary-time derivative. Its integral depends only 
on the initial and final states and hence has no effect on the classical equations of motion, 
but yields a boundary contribution to the Euclidean action. However, it was shown that 
this term, known as the topological phase term, is of central importance for the quantum 
interference effect and makes the tunneling behaviors of integer and half-integer excess spins 
strikingly different 

The Euclidean transition amplitude from an initial state \9i, (pi) to a final state \6f, < 
can be expressed as the following imaginary-time path integral in the spin-coherent-state 
representation, 

(6 f , <Pf\e- nT \O u = J V{0}V{<P} ex P [-S E (6, <f>)}, (3) 

where the Euclidean action Se{@,4>) has been defined in Eq. (2). In the semiclassical limit, 
the dominant contribution to the transition amplitude comes from finite action solutions of 
the classical equations of motion (instantons). According to the standard instanton tech- 
nique, the tunneling rate F for MQT or the tunnel splitting A for MQC is given by T(or 
A) = y4e _,Sci ,il where S c i is the WKB exponent or the classical action which minimizes the 
Euclidean action of Eq. (2). The preexponential factor A originates from the quantum 
fluctuations about the classical path, which can be evaluated by expanding the Euclidean 
action to the second order in small fluctuations.^ It is noted that the above result is based 
on tunneling at the ground state, and the temperature dependence of the tunneling fre- 
quency (i.e., tunneling at excited states) is not taken into account. The instanton technique 
is suitable only for the evaluation of the tunneling rate at the vacuum level, since the usual 
(vacuum) instantons satisfy the vacuum boundary conditions. Different types of pseudopar- 
ticle configurations were developed which satisfy periodic boundary condition (i.e., periodic 
instantons or nonvacuum instantons) .0 

For a particle moving in a double- well-like potential U (x) , the WKB method gives the 
tunnel splitting of nth excited states as@ 

A £ n = ^^exp [-£(£„)], (4) 
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with the imaginary-time action is 



x 2 (E n ) 



S (E n ) = V2m / dxJU (x) - E n , (5) 

where X\^ {E n ) are the turning points for the particle oscillating in the inverted potential 
—U (x). to (E n ) = 2n/t (E n ) is the energy-dependent frequency, and t (E n ) is the period of 
the real-time oscillation in the potential well, 

t(E n ) = V2^ 6 
Jx 3 (E n ) ^E n -U{x) 

where 23,4 (E n ) are the classical turning points for the particle oscillating inside U (x). The 
functional-integral and the WKB method showed that for the potentials parabolic near the 
bottom the result Eq. (4) should be multiplied by 1 2 "e n n! — This factor is very close 
to 1 for all n: 1.075 for n = 0, 1.028 for n = 1, 1.017 for n = 2, etc. Stirling's formula for 
n\ shows that this factor trends to 1 as n — >• 00. Therefore, this correction factor, however, 
does not change much in front of the exponentially small action term in Eq. (4). Recently, 
the crossover from quantum to classical behavior and the associated phase transition have 
been investigated extensively in nanospin systemsSilH^ 

III. MQC FOR BIAXIAL SYMMETRY 

In this section, we consider an AFM system with biaxial symmetry, i.e., which has two 
degenerate easy directions in the biaxial plane. Now the magnetocrystalline anisotropy 
energy is 

E (6, (j)) = K x cos 2 9 + K 2 sin 2 9 sin 2 + E , (7) 

where K\ and K2 are the transverse and longitudinal anisotropic constants satisfying K\ ^> 
K 2 > 0, and Eq is a constant which makes E (6, 0) zero at the initial state. As Ki^> K 2 > 0, 
the Neel vector is forced to lie in the 9 = 7r/2 plane, so the fluctuations of 9 about 7r/2 are 
small. Introducing 9 = n/2 + a (\a\ 1), Eq. (7) reduces to 
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E(a,<j>) » K ia 2 + K 2 sm 2 (f). (8) 

The ground state corresponds to the Neel vector pointing in one of the two degenerate easy 
directions: 9 = n/2, and = 0, it, other energy minima repeat the two states with period 
2n. Performing the Gaussian integration over a, we can map the spin system onto a particle 
moving problem in one- dimensional potential well. Now the transition amplitude becomes 
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with S to t = 2S — s being the total spins of two sublattices, M = jr + 



hS 2 
JV 



being the effective mass, and U (0) = (K 2 V/h) sin 2 being the ef- 
fective potential. J is the exchange density between two sublattices, which is related to 
the transverse susceptibility x± by taking the simple estimate as x± « h 2 j 2 S 2 /JV 2 ^ 
S = miV/hj is the total spin in mi sublattice, and s = mV/hj is the excess spin due to 
the noncompensation of two sublattices. It is noted that the total derivative in Eq. (2), 
when integrated, gives an additional phase factor to the transition amplitude Eq. (9) which 
depends on the initial and final values of 0. For the biaxial symmetry, this phase factor in 
Eq. (9) is exp (—in Stot)- The potential U (0) is periodic with period ir, and there are two 
minima in the entire region 2ir. We may look at U (0) as a superlattice with lattice constant 
7r and total length 2ir, and we can derive the energy spectrum by applying the Bloch the- 
orem and the tight-binding approximation. The translational symmetry is ensured by the 
possibility of successive 2tt extensions. 

The periodic instanton configuration P which minimizes the Euclidean action in Eq. (9) 
satisfies the equation of motion 

--, V 

2 V dr 



)M ^) -U(cf> p ) = -E, (10) 



where E > is a constant of integration, which can be viewed as the classical energy of the 
pseudoparticle configuration. Then we obtain the kink-solution as 



sin (j) p = 1 — k sn (a>iT, fc) , 
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where sn(u;iT, k) is the Jacobian elliptic sine function of modulus k, 
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with oji = V%j[s I i / K j J \ 21 anc ^ n i = \Jk 2 V/?iE > 1. In the low energy limit, i.e., 
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E — > 0, k — > 1, sn(w, 1) — > tanhw, we have 



1 



sin i 



(13) 



p cosh (uj\t) ' 

which is exactly the vacuum instanton solution derived in Ref. 10. 

The Euclidean action of the periodic instanton configuration Eq. (11) over the domain 
(—(3,(3) is found to be 
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where K (k) and E (k) are the complete elliptic integral of the first and second kind, respec- 
tively. Now we discuss the low energy limit where E is much less than the barrier height. 
In this case, k' 2 = 1 — k 2 = hE/K 2 V C 1, so we can perform the expansions of K (k) and 
E (k) in Eq. (15) to include terms like k' 2 and k' 2 In (4//c'), 
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With the help of small oscillator approximation for energy near the bottom of the potential 
well, E = £% a = (n + 1/2) oui, Eq. (15) is expanded as 
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Then the general formula Eq. (4) gives the low-lying energy shift of nth excited states for 
AFM particles with biaxial crystal symmetry at zero magnetic field as 
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When n = 0, the energy shift of the ground state is 
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Now we discuss briefly the dissipation effect and the temperature dependence of the 
decay rate. It is noted that Eqs. (20) and (21) are obtained under the condition that 
the levels in the two wells are degenerate. In more general cases, the transition amplitude 
between two levels separated by the barrier or the decay rate should be sensitive to this 
resonance condition for the two levels. For a spin tunneling problem, it is important to 
consider the discrete level structure. It was quantitatively shown that the phenomenon of 
MQC depends curcially on the width of the excited levels in the right well.Q Including the 
effects of dissipation, the decay rate, in particular, is given bySil 
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where AS n is the level splitting, n' are the levels in the other well and Q nn i is the sum 
of the linewidths of the nth and n'th levels caused by the couplying of the system to the 
environment. For the exact resonance conditions, the temperature dependence of the decay 
rate is 

r(T)=£^^exp(-£: n /3), (23) 

where the level broadening Q n contains all the details of the coupling between the magnet 
and its environment. If the width caused by the coupling of the system to the environment 
is sufficiently large, the levels overlap, so that the problem is more or less equivalent to the 
tunneling into the structureless continuum.0 In this case, the results obtained in this paper 
should be changed by including the dissipation. It is noted that the purpose of this paper is 
to study the coherently quantum tunneling and spin-phase interference at excited levels for 
AFM particles in the absence of a magnetic field at sufficiently low temperatures. Strong 
dissipation is hardly the case for magnetic systems,^ and thereby our results are expected 
to hold. It has been argued that the decay rate should oscillate on the applied magnetic 
field depending on the relative magnitude between the width and the level spacing.i'@0 
However, it is not clear, to our knowedge, what should be the effect of finite temperature in 
the problem of spin tunneling. The full analysis of spin tunneling onto the precession levels 
remains an open problem. 

Now we consider the transition exponent which is usually addressed by experiments. 
Transitions between two states in a bistable system or escaping from a metastable state can 
occur either due to the quantum tunneling or via the classical thermal activation. In the 
limit of T — > 0, the transitions are purely quantum- mechanical and the rate goes as T ~ 
exp (—S c i), with S c i being the classical action or the WKB exponent which is independent of 
temperature. As the temperature increases from zero, thermal effects enter in the quantum 
tunneling process. If the temperature is sufficiently high, the decay from a metastable 
state is determined by processes of thermal activation, and the transition rate follows the 
Arrhenius law, T ~ exp (— U/ksT), with ks being the Boltzmann constant and U being the 
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height of energy barrier between the two states. Because of the exponential dependence of 
the thermal rate on T, the temperature T c characterizing the crossover from quantum to 
thermal regime can be estimated as /csT c = U/Sd- For the present case, one can estimate 
that 



k B T c = 




characterizing the crossover from quantum to classical regime. Typical values of parameters 
for single-domain AFM nanoparticles are: the sublattice spin S = 5000, the excess spin 
s = 150, the longitudinal anisotropy constant K2 ~ 10 5 erg/cm 3 , the transverse anisotropy 
constant K\ ~ 10 6 erg/cm 3 , the exchange energy density between two sublattices J ~ 10 9 
erg/cm 3 , and the radius of particle is about 5 nm. By taking these values, we obtain that 
T c ^704 mK, which agrees well with the experimental result for the diluted samples of 
horse-spleen ferritin J! 

It is noted that hA£^ m is only the level shift induced by tunneling between degenerate 
excited states through a single barrier. The effective periodic potential U (0) = U (<f) + nn) 
can be regarded as a one-dimensional superlattice with lattice constant ir. The tunneling 
through one barrier leads to the level splitting which extends formally to an energy "band" by 
translation symmetry ( the rotation symmetry in the present case, U (<f)) = U (<p + nir)). The 
energy "band" structure of this problem is formally the same as that of a one-dimensional 
tight-binding model in solid state physics. Then the energy spectrum of low-lying excited 
levels can be determined by the Bloch theorem. It is easy to show that if S^ a are the 
degenerate eigenvalues of the system with infinitely high barrier, the energy level spectrum 
is given by the following formula with the help of tight-binding approximation, 

E Ua = £ Ua _ 2Ag Ua cog ^ + Q ^ (24) 

The Bloch wave vector £ can be assumed to take either of the two values and 1 in the first 
Brillouin zone. It is noted that in Eq. (24) we have included the contribution of topological 
phase for AFM particles with biaxial crystal symmetry (i.e., nS to t)- One can easily show 
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that the low-lying tunneling level spectrum, which corresponds to the splittings of nth 
excited state due to the resonant quantum coherence of the Neel vector between energetically 
degenerate states, depends on the parity of S to t (or the excess spin s) significantly. 

At the end of this section, we discuss the possible relevance to the experimental test 
for spin-parity effects in single-domain AFM nanoparticles. One can easily show that the 
specific heat for integer excess spins is much different from that for half-integer excess spins 
at sufficiently low temperatures T ~ T = hAE^ 0, jk B . When the temperature is higher 
hASo <C k B T < Hcui, the excited energy levels may give contribution to the partition 
function of tunneling states. Now the partition function without the dissipation is given by 



1 + - (l — e-^ 1 ) ((3hAe b ta ) 2 1 (2qe'^ l/2 ) 



(25) 



2 

for both integer and half-integer excess spins. Zq = 2e _/37kJl / 2 / (l — e~P hu)1 ^j is the partition 
function in the well calculated for k B T <C AU over the low-lying oscillator like states with 
gbia _ ^ n _|_ i/2) c^!. I (x) = J2 n =o{ x /2) 2n / ( n ') 2 i s the modified Bessel function, and q is 
shown in Eq. (21). We define a characteristic temperature T that is solution of equation 
q e -hw!/2k B T _ i Then we obtain the specific heat up to the order of ((3hA£ ) 2 as 

c = k B (Phu,) 2 — _ + -k B ((3hAS ) 2 { [2 (l - e"^ 1 ) + 4 {(3%^) e~^ 

- (Phuj!) 2 e-^Io (2qe-^ 2 )} - q (ptuv,) [± (5e~ 3 ^/ 2 - e'^ 2 ) 

+4 _ j/ (2 qi e-^' 2 ) + q 2 (Phcu,) 2 (eT^I 2 - e" 3 ^/ 2 ) 

x/^(2ge- /3 ^ l/2 )}, (26) 

for both integer and half-integer excess spins, where I' Q = — ii, and Iq = I 2 — h/x. l v (x) = 
J2 n =o ( — I)™ {x/2) 2n+u /n\T (n + v + 1), where T is Gamma function. The results show that 
the spin-parity effect will be lost at high temperatures. The specific heat for integer excess 
spins is almost the same as that for half-integer excess spins. 
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IV. MQC FOR TRIGONAL, TETRAGONAL AND HEXAGONAL SYMMETRIES 



In this section, we will apply the method in Sec. Ill to study resonant quantum tun- 
neling of the Neel vector in AFM particles with trigonal, tetragonal and hexagonal crystal 
symmetry. For the trigonal symmetry, 



E (9, 0) = K x cos 2 9 - K 2 sin 3 9 cos (30) + E , 



(27) 



where K\ ^> K 2 > 0. The energy minima of this system are at 9 = ir/2, and = 0, 27r/3, 
47r/3, and other energy minima repeat the three states with period 2ir. The spin tunneling 
problem can be mapped onto a problem of one- dimensional motion by integrating out the 
small fluctuations of 9 about tt/2, and for this case U (0) = 2 (K 2 V/h) sin 2 (30/2). Now 
U (0) is periodic with period 27r/3, and there are three minima in the entire region 2tt. The 
periodic instanton configuration with an energy E > is sin 2 (§0 P ) = 1 — k 2 sn 2 (uj 2 T,k), 

r, and m = ^2K 2 V/hE > 1. The low 



where k 



1) /n 2 , u 2 = 3^ 



K 2 J 



1+ 
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energy limit of this periodic instanton configuration agrees well with the vacuum instanton 
solution obtained in Ref. 18. The associated classical action is S p = W + 2E(3, with 



E{k) - (l -k 2 )K(k) 



(28) 



The general formula Eq. (4) gives the low-lying energy shift of nth excited state as 

n+l/2 



HAS. 



tri 



32 



2irn\ 



K2 



i + 
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x exp 
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1 + 



J 



(29) 



2 VXi, 

The periodic potential U (0) can be viewed as a superlattice with lattice constant 27r/3 and 
total length 2n, and the Bloch theorem then gives the energy level spectrum of nth excited 
state £*" = {n+ 1/2) u 2 as £*" = f*" - 2A£% i cos [(-5 tot + 2tt/3], where ^ = -1,0,1 
in the first Brillouin zone. It is easy to show that the low-lying energy level spectrum is 
h£%" l —2hA£%" 1 , and h£ l " +hA£ t r [ l for integer excess spins, the latter being doubly degenerate. 
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While the level spectrum is US^ 1 — hA£^\ and ft£*" + 2hA£%" 1 for half-integer excess spins, 
the former being doubly degenerate. 

For the tetragonal symmetry, the magnetocrystalline anisotropy energy is 



E (9, 0) = K x cos 2 9 + K 2 sin 4 9 - K' 2 sin 4 9 cos (40) + E , 



(30) 



where K\ 3> K 2 , K' 2 > 0. The energy minima are at 9 — ir/2, and = 0, 7r/2, 7r, 37r/2, and 
other energy minima repeat the four states with period 2-7T. The problem can be mapped 
onto a problem of particle moving in one-dimensional potential U (0) = 2 (K' 2 V/K) sin 2 (20) 
by integrating out the small fluctuations of 9 about n/2. Now U (0) is periodic with period 
7r/2, and there are four minima in the entire region 2n. The periodic instanton configuration 
with an energy E > is sin 2 (20 p ) = 1 — k 2 sn 2 (uj 3 t, A;), where k = \J(n\ — 1) /n 2 , cj 3 = 
: , and ni = ^2K' 2 V /HE > 1. The associated classical action is <S P = 



K' 2 J 



W + 2£/3, with 

W = 2S< 



\ J 



1 + 



J 



-V 

SJ 



2 \K X 

The low-lying energy shift of nth excited state is 

2 3/2 



E (k) - (l - k 2 ) K (k) 
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KUV 



^"W 1+ K*)(f) ; 

x exp —25, 



85, 



n+l/2 



1 + 



J 



(32) 



2 VXi, 

Now C/ (0) can be viewed as a superlattice with lattice constant rc/2 and total length 2tc, and 
the Bloch theorem gives the energy level spectrum of nth excited state £„ = (n + 1/2) lu 3 as 
E l n e = £* e - 2A£* e cos [(Stat + tt/2], where £ = -1, 0, 1, 2 in the first Brillouin zone. Then 
the low-lying energy level spectrum is h£^ ± 2hA£^, and ft£* e for integer excess spins, the 
latter being doubly degenerate. While the level spectrum is %£^ ± y/2hA£^ with doubly 
degenerate for half-integer excess spins. 
For the case of hexagonal symmetry, 



E (9, 0) = K x cos 2 9 + K 2 sin 4 9 + K 3 sin 6 9 - K' 3 sin 6 9 cos (60) + E , 



(33) 
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where K\ 3> K 2 , K 3 , K' 3 > 0. The easy directions are at 9 = n/2, and = 0, 7r/3, 2tc/3, tt, 
47r/3, 57r/3, and other energy minima repeat the six states with period 2n. For the present 
case, U ((f)) = 2 (K' 3 V/K) sin 2 (30) is periodic with period ir/3, and there are six minima 
in the entire region 2n. The periodic instanton configuration at a given energy E > 
is sin 2 (30 p ) = 1 — k 2 sn 2 (u; 4 t, k), where k = \J (n 2 — 1) /n 2 , cj 4 = 6^^/ — ~ ^j^ 3 y an d 

n x = ^2K' 3 V /TiE > 1. Correspondingly, the classical action is 5 P = + 2£'/3, with 

E (fc) - (l-k 2 )K{k)] , 



4, K' 2 
3 \T 



1+ 5(^)(f 



(34) 



and the low-lying energy shift of nth excited state is 



%A£ he = 



x exp 
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S\l 



2 (xj {s 



K> 



1 + 



1 / J 



— v n+1/2 



^ J [ 2 Vi^i 



i 5 
3 



(35) 



Now [/ (0) can be regarded as a one-dimensional superlattice with lattice constant n/3. By 
applying the Bloch theorem and the tight-binding approximation, we obtain the energy level 
spectrum of nth excited state E h a e = (n + 1/2) uj a as E h n e = £% e - 2A£ h a e cos [(Sua + f) tt/3], 
where £ = —2,-1,0,1,2,3. If the excess spin s is an integer, the low-lying energy level 
spectrum is h£% e ±2HA£% e , and h£^ e ±hA£^ e , the latter two levels being doubly degenerate. 
If s is a half-integer, the level spectrum is h£^ e ± y/3hA£^ e , and fr£^ e , all three levels being 
doubly degenerate. 

In brief, the low-lying energy level spectrum for trigonal, tetragonal and hexagonal sym- 
metry are found to depend on the parity of the excess spins of AFM particles distinctly, 
resulting from the phase interference between topologically different tunneling paths. The 
structure of low-lying tunneling level spectrum for the trigonal, tetragonal, or hexagonal 
symmetry is found to be much more complex than that for the biaxial symmetry. The 
low-lying energy level spectrum can be nonzero even if the excess spin is a half-integer for 
the trigonal, tetragonal, or hexagonal symmetry. The results of AFM nanoparticles with 
general structure of magnetocrystalline anisotropy will be helpful for experimental test. 
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V. CONCLUSIONS 



In summary, we have studied the topological phase interference effects in the model for 
mesoscopic AFM particles with uncompensated excess spin for the more general structure 
of magnetic anisotropy, such as biaxial, trigonal, tetragonal, and hexagonal crystal symme- 
tries. The low-lying tunnel splittings between nth degenerate excited states of neighboring 
wells are evaluated with the help of the periodic instanton method, and the energy level 
spectrum is obtained by applying the Bloch theorem and the tight-binding approximation 
in one-dimensional periodic potential. This is the first complete study, to our knowledge, of 
spin-phase interference between excited-level tunneling paths in AFM particles with general 
structure of magnetocrystalline anisotropy, which will be useful for experimental check. 

One important conclusion is that for all the four kinds of crystal symmetries, the low- 
lying energy level spectrum for integer excess spins is significantly different from that for 
half-integer excess spins, resulting from the phase interference between topologically distinct 
tunneling paths. For AFM particles with simple biaxial symmetry, which has two degen- 
erate easy directions in the basal plane (i.e., the double- well system), the tunnel splitting 
is suppressed to zero for half-integer excess spins due to the destructive phase interference 
between topologically different tunneling paths connecting the same initial and final states. 
However, the structure of low-lying tunneling level spectrum for the trigonal, tetragonal, 
or hexagonal symmetry is found to be much more complex than that for the biaxial sym- 
metry. The low-lying energy level spectrum can be nonzero even if the excess spin is a 
half-integer for the trigonal, tetragonal, or hexagonal symmetry. Our analytical study pro- 
vides a nontrivial generalization of Kramers degeneracy for double- well system to coherently 
spin tunneling at ground states as well as low-lying excited states for AFM systems with 
m-fold rotational symmetry around the z axis. Note that these spin-parity effects are of 
topological origin, and therefore are independent of the magnitude of excess spins of AFM 
particles, the shape of the soliton and the tunneling potential. One can easily show that the 
heat capacity of low-lying magnetic tunneling states depends significantly on the parity of 
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excess spins for AFM particles with different symmetries at sufficiently low temperatures, 
providing a possible experimental method to examine the theoretical results on topological 
phase interference effects. Our results presented here should be useful for a quantitative un- 
derstanding on the topological phase interference or spin-parity effects in resonant quantum 
tunneling in single-domain AFM particles with different symmetries. 

Over the past years a lot of experimental and theoretical works were performed on 
the spin tunneling in molecular M1112- Ac§ and FegS clusters having a collective spin state 
S = 10 (in this paper S = 10 3 -10 4 ). More recently, Wernsdorfer and Sessoli@ have measured 
the tunnel splittings in the molecular Fes clusters, and have found a clear oscillation of the 
tunnel splitting with the field along hard axis, which is a direct evidence of the role of the 
Berry phase in the spin dynamics of these molecules. Further experiments should focus on 
the level quantization of collective spin states of S = 10 2 -10 4 and their quantum spin phases. 
The theoretical calculations performed in this paper can be extended to the FM and AFM 
particles in a magnetic field. Similar spin-phase interference effects observed in Fe§ cluster 
should be found in single-domain AFM nanoparticles in a magnetic field. Work along this 
line is still in progress. With current technology and fast progress on this field, our study 
on spin-phase interference and resonant quantum coherence effects in AFM nanoparticles 
should be experimentally testable in the near future. 
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